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Abstract. Recently we have proposed a new cut-off scheme for pion loop integrals in the two-pion exchange 
potential. This method allows for a consistent implementation of constraints from pion-nucleon scattering 
and has been successfully applied to peripheral nucleon-nucleon partial waves. We now consider low partial 
waves in the non-perturbative regime, where the regularized Lippmann-Schwinger equation has to be 
solved in order to generate the bound and scattering states. We observe an improved description of most of 
the phase shifts when going from next-to- to next-to-next-to-leading order in the chiral expansion. We also 
find a good description of the deuteron properties. In addition, the new cut-off scheme allows to avoid the 
presence of unphysical deeply bound states. We discuss the cut-off dependence of the four-nucleon low- 
energy constants and show that their numerical values can be understood in terms of resonance saturation. 
This connects the effective field theory approach to boson exchange phenomenology. 

PACS. 13.75.Cs Nucleon-nucleon interactions - 21.30.-x Nuclear forces - 12.39.Fe Chiral Lagrangians 



1 Introduction 

, In Ref. £Q (which is called I form here on), we have pre- 
sented a method to improve the convergence of the chiral 
expansion for the nucleon-nucleon (NN) interaction based 
on spectral function regularization. In earlier approaches, 
an unphysically strong attraction in the isoscalar central 
part of the chiral two-pion-exchange (TPE) at next-to- 
next-to-leading (NNLO) order in the chiral expansion was 
found. This is due to the high-momentum components 
of the exchanged pions, which appear when using dimen- 
sional regularization (or equivalent schemes), and which 
cannot be properly treated in the corresponding effective 
field theory (EFT). Using a cut-off (or spectral function) 
regularization instead of the dimensional one and taking 
reasonable values for the momentum space cut-off allows 
to remove spurious short-distance physics associated with 
high-momentum intermediate states and to greatly im- 
prove the convergence of the chiral expansion. In partic- 
ular, one can use without problems the values of the di- 
mension two low-energy constants (LECs) Ci consistent 
with elastic pion-nucleon scattering data. More precisely, 
in I we have considered the spectral functions obtained 
from the next-to-leading order (NLO) and NNLO TPE 
contributions and argued that only masses below the chi- 
ral symmetry breaking scale should be taken into account 
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explicitly in the loop integrals while shorter range con- 
tributions have to be represented by contact interactions. 
This can be easily implemented by applying a cut-off to 
the spectral functions. We have also proposed a simple and 
convenient way to derive analytic expressions for regular- 
ized TPE in the momentum space based on the spectral 
function representation. 



In I, we have considered the peripheral partial waves 
(I > 2), because at NNLO, these are given entirely by one- 
pion-exchange (OPE) and TPE with no free parameters. 
We have calculated these phases in Born approximation 
which should be legitimate at least for the D- and higher 
waves. The results for the D- and F-waves are still not 
completely converged at NNLO, but the error of a few (1)° 
at i^iab = 300 MeV for the D- (F-)waves appears rea- 
sonable. There is no breakdown of the chiral expansion 
for D-waves beyond Ti a b = 50 MeV and for F-waves be- 
yond Ti a b = 150 MeV as found earlier using dimensional 
regularization. In this paper, we apply our method (cut- 
off regularization, CR, for short) to the low partial waves 
in the non-perturbative regime, where we have to solve 
the regularized Lippmann-Schwinger equation to generate 
the bound and scattering states. As will be demonstrated, 
there are no deeply bound states (for a reasonable range 
of cut-offs in the spectral function representation of the 
effective potential and regularized Lippmann-Schwinger 
equation) , and low-energy observables are not affected by 
CR. 
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Our manuscript is organized as follows. In section 
we briefly review the formalism detailed in I and discuss 
the regularization of the Lippmann-Schwinger equation. 
In section OH we apply the formalism to np partial waves, 
with particular emphasis on the low phases (S- and P- 
waves). We also discuss the range parameters in the S- 
waves and the deuteron properties and compare to the 
results obtained in dimensional regularization (DR). The 
physics behind the values of the low-energy constants ac- 
companying the four-nucleon operators is also elucidated. 
The summary and conclusions are given in section 0] 
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as well as short-distance contact interactions (local 
fermion terms with two derivatives) 
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In I we have calculated the NN potential at NNLO in the 
low-momentum expansion using cut-off regularization in 
the spectral function representation. For completeness, we 
briefly review the pertinent formalism here. 

The chiral TPE potential up to NNLO can be decom- 
posed into isoscalar and isovector central, spin-spin and 
tensor components, generically called Wi(g) here. These 
functions Wi(q) can be represented by a continuous su- 
perposition of Yukawa functions (modulo subtractions) 



Wi(q) = - 



dfj, n 



2M» 



(2.1) 



where the <ii (/z) are the corresponding mass spectra (spec- 
tral functions). Further, q is the momentum transfer in the 
centre-of-mass system (c.m.s.), i.e. q = p' — p, where p' 
and p are final and initial nucleon momenta, respectively, 
and q = |q|. These spectral functions contain the whole 
dynamics related to the exchanged two-pion system and 
can be obtained from the potential via |2] 



<n(n) =Im W;(0 + -in) 



(2.2) 



The spectral function regularization proposed in I sup- 
presses the large-/z contributions to the integrals Eq. (|2.1H 
via a sharp cut-off 



(2.3) 



and thus regulates the short-distance contributions of the 
TPE in a natural way. 

Next, we briefly review the chiral expansion of the NN 
potential. The LO potential is given by OPE and two 
contact interactions, — V}®' + V^t> w i tn 
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where a (t) are nucleon spin (isospin) matrices. The NLO 
corrections are due to two-pion exchange 
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Here, k = (p' + p)/2, and the cut-off regularized loop 
function L A (q) reads: 



L A {q) = 6(A-2M 7r ) — In 
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The regularized expression for TPE provides an explicit 
exclusion of the short-range components in the spectrum 
(i.e. those ones with the range r < A^ 1 ). Furthermore, at 
NLO one also has a correction to OPE, It takes the form 



VP 
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where the LEC dx$ is related to the Goldberger-Treiman 
discrepancy. Thus, the complete NLO contribution is given 



by 2 



(2) 



Finally, the cut-off regular- 



ized NNLO corrections are represented by the subleading 
TPE potential. It takes the form: 
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In what follows, we use these values for the pion de- 
cay constant F n , the pion masses M„± , M n o and the nu- 
cleon mass to: F„ = 92.4 MeV, M n ± = 139.570 MeV, 
M T o = 134.977 MeV, to = 938.918 MeV. In I we adopted 
the following values for the nucleon axial-vector coupling 
g A and the LEC d 18 : g A = 1.26, d ls = -0.97 GeV~ 2 . Al- 
ternatively, one can use the larger value gA = 1.29 and 
completely neglect the NLO correction to OPE given in 
Eq. (|2.8|l (i.e. set dig = 0). In this work we will adopt 
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this second possibility. Notice, however, that such a re- 
placement is not valid in a general case, since the corre- 
sponding chiral gA~ and ciis-vertices with three and more 
pion fields are different. For the LECs ci^ we adopt the 
central values from the Q 3 -analysis of the irN system 
0: ci = -0.81 GeV-\ c 4 = 3.40 GeV" 1 . For the con- 
stant C3 the value C3 = —3.40 GeV -1 is used, which is 
on the lower side but still consistent with the results from 
Ref. 0: c 3 = -4.69 ± 1.34 GeV" 1 . This value of the LEC 
C3 was found in Ref. 0j to be consistent with empirical 
NN phase shifts as well as the results from dispersion and 
conventional meson theories. Further, the same values of 
the LECs 01,3.4 will be used in the upcoming N 3 LO anal- 
ysis which will be published separately [S]- Interestingly, 
similar values for the LEC c 3 have been extracted recently 
from matching the chiral expansion of the nucleon mass 
to lattice gauge theory results at pion masses between 
500 and 800 MeV, see [H] for more details. Notice that at 
NNLO in the chiral expansion of the NN potential one 
could, in principle, use the values of the LECs obtained in 
the Q 2 -analysis of ttN scattering. This certain freedom in 
choosing the values of Ci results in some uncertainty in ob- 
servables which might be viewed as an estimation of some 
higher order effects. We would also like to remark that a 
new np and pp partial wave analysis of the Nijmegen group 
|Zj leads to c 3 = -4.78 ±0.10 GeV^ 1 and c 4 = 3.96 ±0.22 
GeV" 1 using c\ = —0.76 GeV" 1 as input. These values of 
the LECs 03,4 are close to the ones of Ref. [3]. 

Using this potential, one can now generate bound and 
scattering states. For that, consider the partial-wave pro- 
jected Lippmann-Schwinger (LS) equation for the NN T- 
matrix: 

K'Ap'. p) = V«(s>\ P) + T,J §$ V LrlJ>', P") 

V-(^ + i e T '">"'' ) ' (2 ' n) 

where V = + + , m is the nucleon mass. The 
on-shell S- and T-matrices are related via 

S#,(p, p) =S lv - -^pmTfrip, p) . (2.12) 

Contrary to our previous work I, where we have been in- 
terested only in peripheral NN scattering and thus calcu- 
lated the T-matrix perturbatively (i.e. keeping only the 
Born term in Eq. H2.11|l L we now have to solve the LS 
equation non-perturbatively. 

Although we have regularized the TPE contributions 
by cutting off the large-mass components in the spectrum 
(or, equivalently, by explicitly shifting the corresponding 
short-distance components to contact terms), the result- 
ing potential still behaves incorrectly at large q (or equiva- 
lently at small r) . The effective potential is valid for small 
values of the momentum transfer q and becomes mean- 
ingless for q > A x . Moreover, since the potential V grows 
with increasing momenta q, the LS equation (|2.11(l is ul- 
traviolet divergent and needs to be regularized. Follow- 
ing the standard procedure, see e.g. jHJ, we introduce an 



additional cut-off in the LS equation by multiplying the 
potential V(p, p') with a regulator function f A , 

V( P , P')^ f A ( P )V(p, p')f A ( P '). (2.13) 

In what follows, we use the exponential regulator function 

/V) = cxp[-/M 6 ] . (2.14) 

Certainly, both cut-offs A and A are introduced in order 
to remove high-momentum components of the interact- 
ing nucleon and pion fields. The physical meaning and 
the implementation of the cut-offs A and A is, however, 
quite different from each other: while the first one removes 
the short-distance portion of the TPE nuclear force, the 
second one guarantees that the high momentum nucleon 
states do not contribute to the scattering process. One ad- 
vantage of the method proposed in I is that one can (but 
does not have to) choose similar procedures for regulating 
the spectral functions and the LS equation. For further 
discussion of the role and optimal choice of the cut-off A 
in the LS equation the reader is referred to Refs. [9"lll0| . 

In what follows we will vary the cut-off A in the LS 
equation in the range 450 .. . 600 MeV at NLO and 450 .. . 
650 MeV at NNLO, which is a significantly larger range 
than in Ref. The cut-off A is varied independently 
in the range 500 . . . 700 MeV, which is consistent with 
the variation of A. Notice that in principle, more ele- 
gant regularization prescriptions, like e.g. lattice regular- 
ization, would allow to regularize pion loop integrals and 
the Lippmann-Schwinger equation in the same way with- 
out introducing two independent scales A and A. 



3 Results 

For any choice of the cut-offs A and A, the LECs Cs,t and 
C1...7 are fixed from a fit to the np S- and P-waves and 
the mixing parameter e\ for laboratory energies below 100 
MeV. This is the same procedure as used in Ref. In 
the following, we display and discuss our predictions for 
the phase shifts at higher energies and higher angular mo- 
menta as well as for various deuteron properties. We also 
discuss the values of the pertinent low-energy constants 
(LECs) and their interpretation in terms of resonance sat- 
uration along the lines of Ref. ^2] ■ 

Before presenting the results we would like to make 
several comments concerning the theoretical uncertainty 
and our way of estimating it. When performing calcula- 
tions within chiral EFT up to a certain order 0(Q n / X n ), 
where Q ~ M n refers to a generic low-momentum scale 
(soft scale) and A to the scale at which new physics ap- 
pears (hard scale) , the theoretical uncertainty results from 
neglecting the higher order terms and is, in general, ex- 
pected to be of the order ~ (Q/A)™ +1 . In particular, one 
expects for the uncertainty of a scattering observable at 
c.m.s. momentum p to be of the order ~ (max[p, M^]/X) n+ 
One should keep in mind that while the hard scale A is 
governed by A x <~ M p in perturbative calculations in the 
7T7T and ttN sectors, where dimensional regularization is 
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Lab. Energy [MeV] Lab. Energy [MeV] 

Fig. 1. S-wave NN phase shifts versus the nucleon laboratory energy. The light (dark) shaded band shows 
the NLO (NNLO) predictions with CR chiral TPE. The cut-off A in the Lippmann-Schwinger equation is 
varied in the range A = 450 . . . 600 MeV at NLO and A = 450 . . . 650 at NNLO. The cut-off A in the pion 
loops is varied independently in the range A = 500 . . . 700 MeV at both NLO and NNLO. The filled circles 
depict the Nijmegen phase shift analysis (PSA) results 



usually applied and no finite momentum space cut-off is 
introduced, A — minL4, A, M p ] = 450 MeV should be 
adopted in our case. 

Cut-off variation became a common practice to esti- 
mate the uncertainty and to check consistency of non- 
perturbative EFT calculations of few-nucleon systems. 
Low-energy observables should not depend on the cut- 
off value if all terms in the EFT expansion are included. 
In practice, however, calculations are performed at a finite 
order, so that some (small) residual dependence of observ- 
ables on the cut-off remains. One, in general, expects that 
this cut-off dependence gets weaker when higher order 
terms are included. Thus, at first sight one expects nar- 
rower bands for scattering observables at NNLO than at 
NLO for the same variation of the cut-off. This, however, 
does not hold true for the following reason: the cut-off de- 
pendence at both NLO and NNLO has to be compensated 
by inclusion of the contact interactions (counter terms) 
of the order C(Q 4 /A 4 ) and higher. The contact interac- 
tions appear only at even orders 0(Q 21 /A 21 ) in the low- 
momentum expansion while pion exchanges contribute, in 
general, at both even and odd orders. Since the same con- 
tact terms enter the expression for the effective potential 
at NLO and NNLO, similar cut-off dependence for ob- 
servables should be expected at these orders. Variation of 
the cut-off does probably not provide an appropriate esti- 
mation of theoretical uncertainty at NLO since it does not 
rely on missing 0(Q 3 /A 3 )-terms, but only on 0(Q i /X 4 )- 
corrections. Notice further that we were only able to vary 
the cut-off A in the LS equation at NLO in the smaller 
range compared to NNLO, which partially explains why 
the NLO bands in many cases even turn out to be nar- 
rower than NNLO ones. Another reason for that behavior 
has already been discussed in I: variation of the spectral 
function cut-off A has only small effect at NLO, since the 
leading TPE in most cases provides a very small correc- 



tion to the LO potential. The corrections from sublcading 
TPE at NNLO are significantly larger in magnitude (at 
large q), which leads to larger variation of the potential 
associated with the spectral function regularization, see I 
for more details. 



3.1 S-waves 

The phase shifts in the and 3 Si partial waves are 
shown in Fig. ^ One observes a clear improvement when 
going from NLO to NNLO and the description of both 
phases is satisfactory up to the considered energy = 
200 MeV. One should keep in mind that this improvement 
is entirely due to inclusion of the subleading TPE poten- 
tial, since no new contact operators arise at NNLO and 
thus the number of adjustable parameters is the same at 
NLO and NNLO. 

As expected, the bands at NLO and NNLO are roughly 
of the same width. In the 1 S'o-channel, the band at NNLO 
is even wider than at NLO, the latter, however, does not 
properly describe the data at larger energies. We remind 
the reader at this point that the bands at NLO under- 
estimate the uncertainty of the theory at this order. It 
is comforting to see that the bands at NLO and NNLO 
overlap and that the Nijmegen values of the S-wave phase 
shifts are reproduced at NNLO within the theoretical un- 
certainty. 

It is also of interest to consider the scattering length 
and effective range parameters. The effective range expan- 
sion in the S-waves takes the form: 

pcotiS) = — + \rp 2 + v 2 p 4 + v 3 p 6 + v 4 p 8 + 0{p w ) , 
a 2 

(3.1) 

where p is the nucleon centre-of-mass momentum, a is the 
scattering length, r is the effective range and the Vi are the 
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NLO, CR 


NNLO, CR 


NNLO, DR 


Nijmegen PSA 


a [fm] 


-23.447. . . - 23.522 


-23.497... - 23.689 


-23.936 


-23.739 


r [fm] 


2. 60... 2. 62 


2. 62... 2. 67 


2.73 


2.68 


v 2 [fm 3 ] 


-0.46... -0.47 


-0.48... -0.52 


-0.46 


-0.48 


^3 [fm 5 ] 


4.3... 4.4 


4.0. . .4.2 


3.8 


4.0 


v 4 [fm 7 ] 


-20.7. . . - 21.0 


-19.9... -20.5 


-19.1 


-20.0 



Table 1. Scattering length and range parameters for the So partial wave using the CR 
NLO and NNLO potential compared to the DR results (with A = 1000 MeV) and to the 
Nijmegen phase shift analysis (PSA). The values 112,3,4 are based on the np Nijm II potential 
and the values of the scattering length and the effective range are from Ref. |15|. 



shape parameters. The coefficients in the effective range 
expansion are governed by the long-distance physics as- 
sociated with exchange of pions and thus serve as a good 
testing ground for the convergence of the chiral expan- 
sion 14.. To be specific, let us consider the 1 So partial 
wave. For each value of the cut-off A in the Lippmann- 
Schwinger equation one has to determine the values of 
the LECs Ci 5 , Ci s , which accompany the (partial-wave 
projected) contact operators without and with two deriva- 
tives, respectively. These LECs can be fixed, for instance, 
from the first two coefficients in the effective range ex- 
pansion, i.e. from a and r, so that predictions for the Vi's 
can be made. Alternatively, one can fix them from a fit to 
the phase shift at low energy and then calculate all coef- 
ficients a, r and Vi. We will adopt this second method in 
what follows. 

In Table^we present our results for the effective range 
coefficients in the 1 S'o channel. Already the NLO results 
are in a reasonable agreement with the data (as given by 
the Nijmegen PSA). At NNLO we find an improved de- 
scription for all effective range coefficients. The scattering 
length and effective range are still not exactly reproduced 
at this order if a global fit to data is performed. The un- 
certainty for the scattering length and effective range at 
NNLO resulting from variation of the cut-offs A and A 
is of the order of 0.2 and 0.05 fm, respectively, which is 
much larger than the deviations from the experimental 
values (~ 0.05 and ~ 0.01 fm). The shape coefficients are 
reproduced at NNLO within the theoretical uncertainty. 

We also show the result based on the DR NNLO poten- 
tial with the cut-off A in the Lippmann-Schwinger equa- 
tion taken as A = 1000 MeV. Notice that there are two 
unphysical deeply bound states in each of the 1 Sq and 
3 Si~ 3 Di channels in that case. The results for the effec- 
tive range coefficients are a bit less precise than the ones 
obtained with CR TPE. 

Our results for the effective range parameters in the 
3 Si channel are shown in Table Similarly to the previ- 
ously considered case, we observe at NNLO an improved 
description for all coefficients. The predicted values of the 



scattering length and effective range are close to the ex- 
perimental ones (within 0.1% and 1%, respectively). In 
the case of the 3 Si channel, the DR results with A = 1000 
MeV are slightly more precise than the CR ones. 



3.2 P-waves 

Our results for the P-waves and the mixing angle e\ are 
shown in Fig. While the 1 Pi, 3 Pi and 3 P 2 phase shifts 
are visibly improved at NNLO compared to the NLO re- 
sults, the NNLO results for the 3 Po phase shift disagree 
with the data for energies higher than Ei a b ^100 MeV, 
where the NLO results are in a better agreement. At the 
moment, we do not have an explanation for this disagree- 
ment with the data in the 3 Pq channel. We have checked 
that this is corrected at N 3 LO, where a new counter term 
appears in that partial wave. As in the case of S-waves, 
the bands at NLO and NNLO are of a similar width. The 
theoretical uncertainty in the 1 Pi channel is probably un- 
derestimated by the variation of the cut-offs A, A. 

In general, our NNLO results for the phase shifts based 
upon CR TPE look similar to the ones of Ref. |S] , where 
dimensional regularization has been used to calculate pion 
loops. In the latter case, the TPE potential shows un- 
physically strong attraction at intermediate and short dis- 
tances, see jl llfTj for more details. Although a reasonably 
good description is possible with DR TPE at NNLO, as 
documented in [Hj , unphysical deeply bound states arise in 
D- and lower partial waves and one has serious problems 
with the convergence of the chiral expansion. In particu- 
lar, changing the value of the cut-off in the Lippmann- 
Schwinger equation clearly leads to a strong variation of 
the D-wave phase shifts, where the potential still turns 
out to be very strong. The problem with the convergence is 
manifest, since there are no counter terms to compensate 
this cut-off dependence at NNLO. Using spectral function 
regularization in the pion loops as discussed in I we are 
now able to describe the data equally well as with the DR 
version and in addition: 
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NLO, CR 


NNLO, CR 


NNLO, DR 


Nijmegen PSA 


a [fm] 


5.429... 5.433 


5. 424... 5. 430 


5.416 


5.420 


r [fm] 


1.710... 1.722 


1.725... 1.735 


1.756 


1.753 


v 2 [fm 3 ] 


0.06... 0.07 


0.04... 0.05 


0.04 


0.04 


v 3 [fm 5 ] 


0.77... 0.81 


0.71... 0.76 


0.67 


0.67 


v 4 [fm 7 ] 


-4.3. . . - 4.4 


-4.1... -4.3 


-4.1 


-4.0 



Table 2. Scattering length and range parameters for the 3 Si partial wave using the CR 
NLO and NNLO potential compared to the DR results and to the Nijmegen PSA |16) . 



— one can use the same values for the cut-off in the 
Lippmann-Schwinger equation in the LO and NLO 
versions, which are slightly below the mass of the p— 
meson, 

— one does not have spurious deeply bound states, 

— one has a convergent expansion. 

Notice that we use here the same values for the LECs 
as in Ref . EJ and somewhat smaller in magnitude value for 
the LEC C3. Further, we do not include the (incomplete) 
set of leading relativistic (1/m) corrections as done in that 
paper. 



3.3 D- and selected higher partial waves 

The results for D-waves have already been analyzed in I 
making use of the Born-approximation. In Fig.[3]we show 
our results for D-waves obtained by solving the LS equa- 
tion. As expected, the results are quite similar to the ones 
found in I. The small differences like e.g. slightly different 
shape of the 3 £>2 and C2 phase shifts arise due to introduc- 
tion of the exponential regulator function in Eq. (|2.13l) , the 
exact solution of the LS equation as well as due to slightly 
different value of the LEC C3 adopted in the present work. 
Similar results for the phase shifts in the present work and 
in I confirm the high accuracy of the Born approximation 
in these channels. 

As in the previously considered channels, the bands at 
NLO and NNLO are of a comparable width (with excep- 
tion of the 1 I?2 partial wave), and the NNLO results are 
in a better agreement with the data. 

We show in Fig. 0] selected higher partial waves, which 
also display a very similar behavior to the one observed 
in I. We remind the reader that a significant disagreement 
with the data in the 3 G 5 channel at both NLO and NNLO 
should not be considered as a problem because of the ex- 
ceptionally small value of the phase shift in this particular 
channel (more than 10 times smaller in magnitude com- 
pared to other G-waves). 



3.4 Deuteron properties 

We now turn to the bound state properties. We stress 
that we do not use the deuteron binding energy as a fit 
parameter as it is frequently done but rather adopt the 
same parameters as obtained in the fit to the low phases. 
In Table |31 we collect the resulting deuteron properties, 
in comparison to the DR results with A — 1000 MeV and 
the CD-Bonn potential. First, we note a clear improve- 
ment when going from NLO to NNLO. In particular the 
predicted binding energy deviates by 1%— 1.5% at NNLO 
to be compared with ~2%-2.5% deviation at NLO. Also 
visibly improved are the root-mean-square matter radius 
rd and the asymptotic S-wave normalization strength As- 
For rd, a surprisingly good agreement with the data is ob- 
served: the NLO (NNLO) prediction deviate from the ex- 
perimental value by less than 0.35% (0.25%). Deviations 
from the data for the asymptotic S-wave normalization 
strength A s of -2% (-1%) at NLO (NNLO) are of the 
expected size: for a typical deuteron observable one ex- 
pects the uncertainty at NLO (NNLO) to be of the order 
~ M3/(450MeV) 3 - 3% (~ Af4/(450MeV) 4 ~ 1%). The 
quadrupole moment Qd is the only deuteron property for 
which our predictions seem to disagree with the observed 
value by somewhat larger amount than expected, namely 
by ~3.8%-4.5% at NLO and ~3.8%-5.2% at NNLO. It 
is, however, well known that Qd is rather sensitive to 
short-range physics, see e.g.^7]. The deuteron quadrupole 
moment also turns out to be underpredicted by ^4% in 
modern potential model calculations. Notice further that 
similar values for Qd where found in Ref. |18j where DR 
was used and the pertinent LECs were fine-tuned to the 
deuteron binding energy. 

The DR NNLO predictions shown in Table H are m a 
better (somewhat worse) agreement with the data for As 
(Qd) and of the same quality for the other observables as 
the CR NNLO results. 

3.5 Low-energy constants and resonance saturation 

We are now in the position to confront the LECs deter- 
mined from chiral effective field theory with the highly 
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successful phenomenological/meson exchange models of 
the nuclear force following the lines of Ref. In this 
work we will consider the Bonn-B |19| and Nijmegen 93 
20 potentials, which are genuine one-boson-exchange 
(OBE) models. In these models the long range part of 
the interaction is given by OPE (including a pion-nucleon 
form factor) whereas shorter distance physics is expressed 
as a sum over heavier mesons exchange contributions: 



(3.1) 



Here some mesons can be linked to real resonances (like 
e.g. the p-meson) or are parameterizations of certain phys- 
ical effects, e.g. the light scalar-isoscalar cr-meson is needed 
to supply the intermediate range attraction (but it is not 
a resonance). The corresponding meson-nucleon vertices 
are given in terms of one (or two) coupling constant (s) 
and corresponding form factor (s), characterized by some 
cut-off scale. These form factors are needed to regular- 
ize the potential at small distances (large momenta) but 
they should not be given a physical interpretation. As 
depicted in Fig[3 for nucleon momentum transfer below 
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the masses of the exchanged mesons, one can interpret 
such exchange diagrams as a sum of local operators with 
increasing number of derivatives (momentum insertions). 
This is explained in detail in Ref. 12 . In that work we 
power expanded the short-range part of different phe- 
nomenological potential models and compared the result- 
ing contact operators with the ones in the EFT approach. 
The latter have to be corrected by adding the correspond- 
ing power expanded TPE contributions, which are not 
present in the phenomenological models. We have then 
demonstrated explicitly that the values of the LECs Cj de- 
termined from various phenomenological OBE models are 



close to the values found in EFT at NLO and NNLO. The 
TPE contributions have been calculated in 12 s using di- 
mensional regularization. In that work we have restricted 
ourselves to the NNLO version with numerically small val- 
ues of the LECs 03,4, C3 = —1.15, c 4 = 1.20 GeV -1 , which 
are not consistent with the irN system. Using the values 
of these LECs obtained from the irN system in the TPE 
potential calculated with dimensional regularization (or 
equivalent schemes) leads, as explained before, to conver- 
gence problems due to the strong attraction in the cen- 
tral part of the potential. In particular, the D-wave phase 
shifts turn out to be strongly cut-off dependent. Further, 
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M 





Fig. 5. Expansion of a meson exchange diagram in terms of 
local four-nucleon operators. The dashed and solid lines de- 
note the meson M — p,a,ui, . . . and the nucleons, respectively. 
The blob and the square denote insertions with zero and two 
derivatives, in order. The ellipses stands for operators with 
more derivatives. 



unphysical deeply bound states arise in low partial waves. 
Obviously, no resonance saturation can be established for 
this NNLO version of the potential, which is strongly non- 
phasc-equivalent to the phenomenological OBE models. 
For example, the LECs Cj in this NNLO version are typ- 
ically several times larger in magnitude and differ very 
much from the values at NLO. 

We will now demonstrate how resonance saturation 
works for the chiral NN forces at NLO and NNLO in- 
troduced above, where the new spectral function regular- 
ization has been used to derive the TPE contributions. 
Differently to Ref. [T2|, all LECs c; are now consistent 
with irN scattering. In Table 01 we compare the values 
of the LECs C, at NLO and NNLO with the ones re- 



sulting from the OBE models as explained before. We re- 
mind the reader that the contribution from chiral TPE 
should be accounted for properly in order to allow for a 
meaningful comparison with the OBE models. To achieve 
that we power expand the chiral TPE at NLO and NNLO 
and identify the corresponding contributions to the LECs, 
which are given analytically in appendix fA] Notice that 
differently to 12; these contributions depend now on the 
spectral function cut-off A. The second and third columns 
in Table 0] show the corresponding numerical results at 
NLO and NNLO for the central value A = 600 MeV, re- 
spectively. The indicated uncertainty refers to the cut- 
off variation A = 500 . . . 700 MeV. The fourth and fifth 
columns contain the values of the LECs Cj at NLO and 
NNLO, where the just discussed contributions from TPE 
have already been added. The numbers are presented for 
our central values of A and A, A — 550 MeV and A = 
600 MeV and the uncertainties refer to variations A — 
500 .. . 700 MeV and A = 450... 600 MeV (A = 450... 650 
MeV) at NLO (NNLO). Notice that the uncertainties due 
to the variation of A and A are very large in several cases: 
for Ciso and Ciso at NLO and for C^po and C3.p1 at 
NNLO. Such a strong variation in the values of the LECs 
arises when the cut-off A in the LS equation becomes 
too large and one leaves the plateau-region for the cor- 
responding LEC Ci(A). This situation is exemplified in 
FigEl where we show the yl-dependence of the LEC C^pi 
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NLO, CR 


NNLO, CR 


NNLO, DR 


CD-Bonn 


Exp. 


E d [MeV] 


-2.171 ... - 2.186 


-2.189. . . - 2.202 


-2.230 


-2.225* 


-2.225 


Qd [fm 2 ] 


0.273 . . . 0.275 


0.271... 0.275 


0.270 


0.270 


0.286 


V 


0.0256... 0.0257 


0.0255 . . . 0.0256 


0.0257 


0.0255 


0.0256 


Td [fm] 


1.973... 1.974 


1.970 . . . 1.972 


1.970 


1.966 


1.967 


As [fm" 1 / 2 ] 


0.868 . . . 0.873 


0.874 . . . 0.879 


0.886 


0.885 


0.885 


Pd[%] 


3.46... 4.29 


3.53... 4.93 


6.71 


4.83 





Table 3. Deuteron properties derived from the CR chiral potential at NLO and NNLO compared to 
the DR NNLO results of j>J, one "realistic" potential and the data. Here, Ed is the binding energy, Qd 
the quadrupole moment, r\ the asymptotic D/S ratio, r<j the root-mean-square matter radius, As the 
strength of the asymptotic S-wave normalization and Po the D-state probability. * denotes an input 
quantity. 



LEC 


TPE (NLO) 


TPE (NNLO) 


d (NNLO) 


d (NNLO) 


Bonn B 


Nijm-93 


Ciso 

Ciso 

C3S1 

Casi 

Co. 

Cipi 

C3P0 

C3P1 

Czpi 


-0 04 + 000 

-0 570 +°' 036 
u -°' u -0.022 

0-013™ 
0-638 ™ 
-0 190 + 0012 

_ 67 +0.007 

-0 425 +0 ' 025 
u.^zo _ 014 

Q 246 +0.009 

u.z'io _ 01B 
-0 22 + 000 

u ' uzz -0.000 


-0 04 + 000 
" ,um -0.001 

-0 4 43 + - 078 

-0 04 + - 000 
" ,um -0.001 

-0 4 43 + - 078 

205 +0 - 024 

-0 090+ 0013 
u.uau _ o oog 

nna +0.003 
O- 006 -0.003 

2 47 +0.032 

n , e 1 +0.020 
u - i0i -0.028 


-0 H7+ 2 - 271 

' -0.042 

1 294 +2 ' 873 

-0.322 

fl 1QC +0.025 
U.-LOO _ .o21 

231 + 112 

U.ZOl „o.007 
n n 9 r +0.000 

u.ozo _ . 036 

146 + 005 

923 + 142 
u.yzo _ 103 

-0 260 +0 003 

U.ZDU _o.oo5 

-0 262 +0 032 

U - ZDZ -0.073 


-0 158+ ' 178 

, 91 o +0.408 
l.ZIO -o 084 

-0 137+0017 

u - 10 ' -0.027 

523 +0 197 

U.OZO _ 039 

-0 395+0 00 7 

U.OJO _ .072 

126 +0 023 

U - 1ZD -0.017 

0-920 

-0.108 ™ 
-0.421 t° Zt 


-0.117 
1.276 

-0.101 
0.660 

-0.410 
0.454 
0.921 

-0.075 

-0.396 


-0.061 
1.426 

-0.014 
0.940 

-0.343 
0.119 
0.802 

-0.197 

-0.467 



Table 4. The LECs d at NLO and NNLO compared with the results from the Bonn B 
and Nijmegen 93 OBE potential models. Also shown are contributions from chiral TPE as 
explained in text. The C t are in 10 4 GeV~ 2 and the d in 10 4 GeV" 4 . 



at NNLO for A = 600 MeV. The yl-dependence of C 3P i is 
similar to the one of the three-nucleon force observed in 
Ref. EH- 1 The first branch in Fig. (for A < 730 MeV) 
corresponds to the case of no deeply bound states. For 
larger values of A unphysical deeply bound states arise. In 
that case the situation is similar to the DR NNLO analy- 
sis of Ref. |S] . Clearly, it only makes sense to discuss res- 
onance saturation of the Ci 's in the plateau-region of the 

1 The long-range part of the potential in that reference be- 
haves like 1/r 2 at short distances, while the CR chiral TPE at 
NNLO is even more singular and behaves like 1/r 5 , see I. 



first branch, where they only change modestly and where 
the effective potential is at least not strongly non-phase- 
equivalent to the OBE models. Notice that the strong 
variation of the LECs with the cut-off did not occur in 
Ref. 12 simply because A was varied in a much smaller 
range and never left the plateau-region. 

The last two columns in Table^Jshow the LECs as pre- 
dicted by resonance saturation based upon the Nijmegen 
93 and Bonn B potential models. As in Ref. J21> we ob- 
serve a remarkable agreement between the LEC values 
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Fig. 6. "Running" of the LEC C3P1 with the cut-off A at 
NNLO. The cut-off A in the spectral function representation 
is fixed at our central value A — 600 MeV. 



obtained from fit to NN phase shifts in the EFT approach 
and the ones resulting from the OBE models. 



4 Summary and conclusions 

In this paper we have considered the two-nucleon poten- 
tial in chiral effective field theory, making use of the novel 
method of regularizing the pion loop integrals introduced 
in I. For the low partial waves and the deuteron consid- 
ered here, one has to solve the LS equation with the prop- 
erly regulated potential, cf. Eq. (|2.I3J) . One thus has to 
deal with two different cut-offs, the first one related to 
the spectral function regularization (denoted A) and the 
second one related to the regularization of the potential 
in the LS equation (denoted A). As pointed out, both of 
these cut-offs can be chosen in the same range. We obtain 
the following results: 

1) We have varied the cut-off A in the LS equation in the 
range 450 .. . 600 MeV at NLO and 450 .. . 650 MeV at 
NNLO, which is a significantly larger range than in 
Ref . JT] . The cut-off A is varied independently in the 
range 500 . . . 700 MeV, which is completely consistent 
with the variation of A. 

2) As shown in Fig.^we obtain a modest/satisfactory de- 
scription of the S-waves at NLO/NNLO when fitting 
the Nijmegen PSA for energies below 100 MeV. Note 
that the NNLO result only differs from the NLO one 
by inclusion of subleading two-pion exchange correc- 
tions. The LECs 01,3,4 which enter this subleading TPE 
contributions are consistent with the analysis of pion- 
nucleon scattering in chiral perturbation theory. The 
resulting S-wave scattering lengths, effective ranges 
and higher order range parameters come out in good 
agreement with the ones deduced from the Nijmegen 
PSA, cf. Tables and [5J In contrast to earlier work 
|H] employing dimensional regularization, we have no 
unphysical deeply bound states in the S-waves (and 
any other partial wave). 



3) At NNLO, all P-waves and the mixing parameter e\ 
are well described at NNLO below E Yah ~ 100 MeV 
and improved as compared to their NLO representa- 
tions (one parameter per phase at NLO and NNLO), 
with the exception of the 3 Pq partial wave. We stress 
that at NNNLO this deficiency is cured [5]- 

4) The D- and higher partial waves come out similar to 
the results obtained in Born approximation in I. For 
further discussion, we refer to that paper. 

5) The deuteron properties collected in Table EH are given 
parameter free. We find a slight improved in the quad- 
rupole moment as compared to the calculation based 
on dimensional regularization (as long as the binding 
energy is not used as an input parameter). 

6) The theoretical error bars on the various partial waves 
and the bound state properties come out consistently 
with expectations, see in particular the discussion in 
the beginning of Section 03 

7) We have shown that the numerical values of the four- 
nucleon LECs are in good agreement with the ones 
derived from semi-phenomenological boson exchange 
models once the two-pion-exchange contribution is prop- 
erly accounted for, as detailed in the appendix. This 
strengthens the conclusions of Ref. ^2] and bridges the 
gap between the chiral EFT approaches and more phe- 
nomenological models describing the forces between 
two nucleons. 

It is now of utmost importance to investigate the next 
order, that is the NNNLO corrections, to a achieve a truly 
accurate description of all important partial waves. Work 
along these line is underway [1] (for a first attempt using 
dimensional regularization see |22|L 
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A Reduction of the two-pion exchange 
contributions 

As stated before, we have to add the contribution of the 
TPE to the LECs so as to be able to compare with boson 
exchange potentials. Expanding the explicit expressions 
for the renormalized TPE potential given before in powers 
of q allows for a mapping on the spectroscopic LECs (of 
course, the TPE contains many other contributions, which 
are, however, of no relevance for this discussion). At NLO 
we get 



^yNLO 
°1S0 



1 /^NLO 

3° 3S1 



NLO 
SO 



= 3a 



r NLO 
°3S1 



18 a (1 + 4g 2 A -8g 4 A )M 2 A 2 
(2 + 17g 2 A - 88. 9 l)i 2 

-2(1 + ig 2 A - i 
!)n - (2 + 17g 2 A - 40 5 l)i 2 
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/-IN. 



NLO 



,-rNLO 
°1P1 



riNho 

°3P0 



^NLO 
°3P1 



,-rNLO 
U 3P2 



+ 2(l + 4g 2 A -8g\)M 2 
-MV2ag 4 A A 2 , 
6a[(2 + 17g A -16g A )A 2 

-2(l + 4g 2 A -8g A )M 2 



(A.l) 



-2a 



= 2a 



2a 



(2 + l7g A + 7Ag\)A 2 

-2(l + 4g A -8g 4 A )M 2 
-(2 + 17g 2 A - 61g A )A 2 

+ 2(l + 4g 2 A -8g A )M 2 
-(2 + \7g\ - 7g\)A 2 

+ 2(l + 4 5 2 -8 5 l)Af 2 



where 



a 



Vi 2 - 4M„ 
432 Fi A 3 7T 



(A.2) 



The above expressions coincide exactly with the ones given 
in Ref. |12| in the limit A — > oo. Similarly, we can give the 
additional TPE NNLO contributions to the various LECs: 



/=fNNLO 

^NNLO 
°1S0 



NNLO 



C™ = -36/3(2 Cl - C3 )M 2 ^, 



NNLO 
3S1 



3/3 



/"(IN J 

■»N1 
'1P1 



/-iNNLO 



^NNLO 
U 3P0 



^.NNLO 
°3P1 



^NNLO 
°3P2 



where 



4c 4 A 2 

- 2ci(5i 2 - 2lM w - AM 2 ) 

+ c 3 (lli 2 - 2AM n - AM 2 )] , 
12^2(5 a A 2 , (A.3) 
2/3 - 12c 4 /i 2 + 2ci(5i 2 - 2iM w - 4M 2 ) 

~c 3 (lli 2 -2iAf ;r -4Af 2 )] , 

2ci(5i 2 -2iAf 7r -4M 2 ) 

- c 3 (lli 2 - 2AM n - 4M 2 ) 
2 (3 2 Ci A 2 + 2 Cl (5A 2 - 2AM n - 4M 2 ) 

-c 3 (lli 2 -2iAf 7T -4Af 2 ) 
2ci(5i 2 -2iA/ 7r -4A/ 2 ) 
-c 3 (lli 2 -2ylAf 7r -4M 2 ) 



= 2(3 



8c 4 A 2 



2(3 



2c 4 A 2 



= 



g 2 A M v (A-2M 7r ) 
48 Fi i 3 



(A.4) 



These expressions depend on the dimension two LECs 
c i,3,4 as discussed before and coincide with the ones given 
in Ref. m the limit A — > oo modulo 1/m— corrections, 
which are not considered in the present work. 



Finally, it should be kept in mind that we use here 
the same notation as in our previous work [SJ, according 
to which the contact terms that are nonanalytic in the 
pion mass and result from TPE, see e.g. I, are not shown 
explicitly. In practical applications to the NN system at 
fixed value of Af w such terms cannot be disentangled from 
the zero-range counter terms (contact interactions from 
the Lagrangian) and therefore do not need to be treated 
separately. A more precise treatment is required if the 
pion mass dependence of the nucleon force is studied |23l 

mm 
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